AD-A071  449  GEORGE  WASHINGTON  UNIV  WASHINGTON  0 C PROGRAM  IN  L06--ETC  F/G  15/5 

STATISTICAL  PROBLEMS  IN  the  CONTROL  OF  MULTI-ECHELON  INVENTORY  —ETC(U) 
JUN  79  S ZACKS  N00014-75-C-0729 

UNCLASSIFIED  SERIAL-T-404  NL 


I OF  / 


M1A071449 


r 


THE 

GEORGE 

WASHINGTON 

UNIVERSITY 


O- 

o 

CJ> 


STUDENTS  FACULTY  STUDY  R 
ESEARCH  DEVELOPMENT  FU~f 
URE  CAREER  CREATIVITY  CC 
MMUNITY  LEADERSHIP  TECh 
NOLOGY  FRO 
ENGINEERING  A 
GEORGE  WASH 


D D C 

mm 


JUl  19  19T9 


liyyisira 

D 


n 7 

INSTITUTE  FOR  MANAGEMENT 
SCIENCE  AND  ENGINEERING 


SCHOOL  or  ENGINEERING 
AND  APPLIED  SCIENCE 


THIS  DOCUMENT  MAS  BEEN  APPROVED  FOR  PUBLIC  RELEASE  AND  SALE,  ITS  DISTRIBUTION  IS  UNLIMITED 


MULTI-ECHELON  INVENTORY  SYSTEMS 
by 

S.  Zacks 


Serial  T-404 
15  June  1979 


Tue  George  Washington  University 
School  of  Engineering  and  Applied  Science 
Institute  for  Management  Science  and  Engineering 


Accession  Per  x 

NITS  GKutl  ^>< 
DDC  TAB  C 
Unannounced  L 
Justification 

r 

B' 

I» 

I 

Istrlbutien/ 

Availability  Codas 

Avail  a 
spec: 

nd/ox* 

al 

Program  in  Logistics 

Contract  N00014-75-C-0729 
Project  NR  347  020 
Office  of  Naval  Research 


Ihis  document  has  been  approved  for  public 
sale  and  release;  its  distribution  is  unlimited. 


D D C 


[T^r?r?Df?nn  r?f?r 

7T\ 

W JUL  I9T9 

Ulki^u  urLbl 

MJ 

- D 


NONE 

SECURITY  CLASSIFICATION  OF  THIS  PACE  (HTiAn  Data  En larad) 

REPORT  DOCUMENTATION  PAGE BEFOREDCOMPLETINGNFORM 

I.  REPORT  NUMBER  [2.  GOVT  ACCESSION  NO. I 3.  RECIPIENT'S  CATALOG  NUMBER 


T-404 

♦ title  fmi  a Submit) 


S.  TYPE  OF  REPORT  ft  PERIOO  COVERED 


STATISTICAL  PROBLEMS  IN  THE  QONTROL  OF 
MULTI-ECHELON  INVENTORY  SYSTEMS 


SCIENTIFIC 


6.  PERFORMING  ORG.  REPORT  NUMBER 


ft.  CONTRACT  OR  GRANT  NUMBERfft) 


N00014-7  5-C-07  29 


T2.~  REPORT  OATE 

/ 15  JUNE  1979  / 


11.  CONTROL  .ING  OFFICE  NAME  AND  ADDRESS 

OFFICE  OF  NAVAL  RESEARCH 

CODE  434  v 

ARLINGTON,  VA  22217 


4.  MONITORING  AGFNCY  NAME  ft  AODRESSflf  dillorant  tram  Controlling  Olltca)  IS.  SECURITY  CLASS,  (o I Ihla  report) 


IS «.  DECL  ASSI  FI  CATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thl»  Report) 

DISTRIBUTION  IS  UNLIMITED;  APPROVED  FOR  PUBLIC  SALE  AND  RELEASE. 


17.  DISTRIBUTION  STATEMENT  (of  the  mbatrmct  entered  In  Block  20,  ft  different  from  Report) 


19.  KEY  WORDf  (Continue  on  revoreo  aide  If  ncctiaaiy  and  Identify  by  block  numbar ) 

MULTI-ECHELON  INVENTORY  SENSITIVITY 

BAYES  ADAPTIVE  CONTROL  DEMAND  DISTRIBUTIONS 

ROBUSTNESS 


20.  ABSTRACT  (Contlnum  on  rmvmroo  mlrtm  If  nacaaaary  and  tdantity  by  block  nun  bar) 

— - • Problems  of  the  robustness  of  adaptive  control  procedures  against  deviations 
of  the  empirical  phenomena  from  the  assumed  models  are  of  high  practical  and 
theoretical  importance.  The  present  paper  studies  the  robustness  of  Bayes 
adaptive  control  of  two-echelon  inventory  systems  relative  to  erroneous  assump- 
tions concerning  (i)  the  initial  control  parameters;  (ii)  the  stationarity  of 
the  demand  distribution  and  (iii)  the  nature  of  the  demand  distribution. 


I JAN  73  1473  EDITION  OF  I NOV  Sft  IS  OBSOLETE 
S/N  0103-014-6601  I 


NONE 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Bntorod) 


li 


r 


THE  GEORGE  WASHINGTON  UNIVERSITY 
School  of  Engineering  and  Applied  Science 
Institute  for  Management  Science  and  Engineering 

Program  in  Logistics 


Abstract 

of 

Serial  T-404 
15  June  1979 


STATISTICAL  PROBLEMS  IN  THE  CONTROL  OF 
MULTI-ECHELON  INVENTORY  SYSTEMS 

by 


S.  Zacks* 


Problems  of  the  robustness  of  adaptive  control  procedures  against 
deviations  of  the  empirical  phenomena  from  the  assumed  models  are  of 
high  practical  and  theoretical  importance.  The  present  paper  studies 
the  robustness  of  Bayes  adaptive  control  of  two-ecnelon  inventory  systems 
relative  to  erroneous  assumptions  concerning  (i)  the  initial  control 
parameters;  (ii)  the  stationarity  of  the  demand  distribution  and  (iii) 
the  nature  of  the  demand  distribution. 


Research  Supported  by 

Contract  N00014-75-C-0729 
Project  NR  347  020 
Office  of  Naval  Research 


*Case  Western  Reserve  University 


T-404 


1.  Introductioi. 

In  a series  of  papers  on  two-echelon  multi- station  inventory  systems 
[3, 4, 6, 7]  optimal  ordering  procedures  were  developed  for  systems  modeled 
towards  possible  naval  applications.  In  these  models  we  assumed  that  order- 
ing at  the  lower-echelon  stations  and  at  the  upper- echelon  depot  is  periodi- 
cally done  (nu  the  beginning  of  every  month)  and  that  the  lead  times  are 
fixed,  one  month  from  the  depot  to  the  lower  echelon  stations  and  two  months 
from  the  wholesale  supply  centers  to  the  depot.  We  further  assumed  that  the 
monthly  demand  at  each  lower- echelon  station  is  a random  variable  having  a 
Poisson  distribution,  stationary  in  time.  Moreover,  the  demand  variables  at 
different  lower- echelon  stations  are  mutually  independent.  The  means  of 
the  corresponding  demand  distributions  were  assumed,  however,  to  be  unknown. 
Adaptive  Bayes  ordering  policies  were  developed  in  [4,t>]  and  proved  in  [7] 
to  be  optimal  for  the  specific  cost  functions  and  objectives  specified. 

These  Bayes  procedures  were  based  on  the  assumption  that  the  unknown  means 
of  the  Poisson  distributions,  X^, X^, ...,X^,  are  priorly  independent  having 

prior  gamma  distributions,  «#(— , v.),  with  scale  parameters  t.  and  shape 

Ti 

parameters  v^,  i=l, ...,k.  The  assumption  about  the  prior  independence  of 
X^, ...,X^  was  essential  to  the  development.  On  the  other  hand,  the  assump- 
tion concerning  the  prior  gamma  distribution  of  the  X,  (i=l, ...,k)  is  not 
essential  but  reasonable,  since  the  gamma  distribution  is  a conjugate  prior 
to  the  Poisson  distribution  (see  [1]).  Moreover,  the  gamma  priors  yield 
negative-binoti ial  predictive  distributions  (see  [£]),  which  have  been  found 
empirically  to  fit  well  the  demand  distributions  of  many  classes  of  items  in 
the  Claris  system  (see  [8]).  The  present  study  investigates  certain  robust- 
ness or  sensitivity  questions  related  to  these  Bayes  control  procedures. 
First,  it  is  interesting  to  evaluate  the  sensitivty  of  the 
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system  to  the  choice  01'  the  initial  prior  parameter:  . Hecond,  it  is  important, 
to  investigate  the  effects  on  the  system  of  possibl ' changes  in  the  intensity 
of  the  demand  distributions  (the  V'isson  means)  which  may  occur  at  unknown 
time  points.  Third,  it  is  worthwhile  to  study  the  robustness  of  the  control 
procedures  to  deviations  of  the  demand  distributions  from  the  assumed  i'oisson 


distributions.  Assuming  that  the  demand  dlstribut ions  are  indeed  i’oisson, 

or  very  close  to  Poisson,  we  may  expect  that  the  con* rol  procedures  will  be 

efficient  if  the  choice  of  the  prior  parameters  i 1, ...,k  is 

sueh  that-  the  prior  mean  v.t.  is  close  t.o  tlie  true  value  of  X.  and  the 

11  1 

prior  variance  vT . , is  small . The  'A ayes,  procedures  are  optimal  in  the 
sense  of  minimizing  the  prior  risk  tuider  the  assumed  prior  parameters,  i.ow- 
ever,  for  t.lie  sensitivity  analytic,  we  compare  tin'  procedures  by  simulating 
certain  realisations  of  the  systems  over  a large  numeer  of  periods,  and  com- 
paring the  moving  averages  of  the  actually  scoured  vests.  if  the  demand  pro- 
cess is  stationary  the  sequence  of  moving  averages  is  ergodie,  converging 
to  the  steady -state  expected  monthly  cost..  If  the  demand  process  is  not 
stationary  the  sequence  of  moving  averages  will  fluctuate  and  will  reflect, 
this  lack  of  stability.  However,  wv  will  demonstrate  that  even  in  such  non- 
et.at  ion ary  eases,  the  Hayes  adaptive  control  perform;  better  than  the  optimal 
control  which  assumes  a complete  knowledge  of  the  demand  dist.rlbut.ion.  One 
can  probably  attain  a similar  degree  of  effectiveness  by  non-Hayes  adaptive 
procedures,  like  the  ones  based  on  the  maximum-likelihood  estimation  of  the 
unknown  paramo. or s.  We  have  not,  however,  compared  Hayes  and  non- Hayes 
procedures.  Finally,  for  the  purpose  of  .-.hewing  the  possible  effect,  on  the 


moving-avorjtgos  el'  t.ho  monthly  costs  of  pronounced  deviation  of  the  demand 
distribution  from  the  i’oisson,  wo  studied  the  behavior  of  the  Hayes  control 
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procedure  when  the  demand  distribution  is  a compound  Poisson,  wi th  mean  at 
each  mo n tli  randomly  determined  according  to  a uniform  distribution  on  some 
finite  range  The  general  conclusion  derived  from  the  present  investigation 
is  that  the  Bayes  adaptive  procedures  adjust  quite  fast  to  possible  changes 
in  tiie  underlying  demand  distributions  and  generally  provide  effective  controls. 
We  provided  in  [:t,o]  also  an  approximate  Bayes  procedures.  We  show  that 
these  approximate  procedures  is  effective  only  if  the  initial  state  of  the 
system  is  not  too  unfavorable.  The  reader  is  refered  to  Section  b for 
summary  and  conclusions  based  on  certain  simulations.  In  order  to  reduce 
somewhat  the  amount  of  computations  needed  for  our  investigation  we  have 
restricted  attention  to  one- station  two-echelon  system,  which  will  be 
specif ically  described  in  the  next  section. 


: . The  Iwo- Echelon  System  and  its  Statistical  Control. 

"he  system  under  consideration  is  comprised  of  a lower- echelon  station, 

E,  and  an  upper- echelon  depot,  P.  Customers  arrive  at  random  at  E and 
demand  a random  number  of  units  of  a specified  item.  Let  X ^ ( i=l, 

designate  the  total  number  of  items  demanded  from  K during  the  i-tli  month. 

It  is  assumed  that  X.  has  a Poisson  distribution  with  mean  \.  (i=l, ?, ...), 
and  X^,  X„, ...  are  mutually  independent,  given  ...  . Let  W.  denote 

the  number  of  units  in  stock  at  E at  the  beginning  of  the  i-th  month.  At 
the  beginning  of  the  i-th  month  K issues  an  order  from  P for  \\  items. 
-J.  s Y.  SB.  - 0 . , wliere  B.  denotes  the  total  number  of  units  at  the 
system.  A negative  Y means  that  items  are  sent  back  to  P.  The  lead-time 
for  the  flow  of  stock  between  E and  P is  L^  = 1 month.  At  the  beginning 
of  the  i-th  month  the  depot,  P,  issues  an  order  of  units  to  be  pur- 

chased from  outside  sources.  It  is  assumed  that  the  lead-time  for  this 
replenishment  is  L.,  P months.  The  cost  of  procurement  of  new  luiits  is 


C*[ $/unit  J.  Ail  item  demanded  at  E can  be  supplied  to  the  customer  at  any 
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time  during  the  month.  Unsatisfied  demand  at  the  end  of  a month  is  lost, 
with  penalty  for  shortage  of  p[$/unit]j  (no  backlogging) . Items  which  are 
left  at  E undemanded  accrue  a carrying  cost  of  0[  $/ item] . The  variables 

i*  i-l**'i^  * l«*«*>k,  designate  the  state  of  the  system  during 

the  i-th  month.  The  following  recursive  relations  hold: 


(T.l) 


Vi = u'i + Yi  • V 


l'i+l  ' (;'i  + "i-l  " ^ * 


i = 0,1, 


where  S . and  (h  are  the  initial  stock  levels,  X = 0,  Y,  e 0, 


'-1 


and 


-0  ° 


are  the  prior  pending  procurement  orders,  and  a = max(0,a). 


In  ( U , r ] we  suggested  to  determine  Y.  in  a manner  which  minimises 
the  expected  cost  to  the  lower- echelon,  duo  to  surplus  or  shortages  during 
the  ( i+1 ) st.  month.  Furtliermore,  it  wa s suggested  that  the  ordering  policy 
of  the  upper- echelon  will  be  the  minimal  required  for  guaranteeing  that  the 
probability  oi  shortage  at  the  P will  not  exceed  a specified  risk  level  (3, 

0 < d < 1.  As  shown  in  [ k ],  if  X^  = = X,  then  the  "desired" 

ordering  level  of  E is  Y^(\, 0 ) = kL\\)  - 0,  where 


(2.2) 


k^ (X ) least  non- negative  integer  k 
satisfying:  Pos(k|2\)  - R • Pos(k|\)  s 0 , 


where  R = p/ ( .'+p)  and  l\'s(k|n)  designates  the  c.d.f.  of  a Poisson  distri- 
bution with  mean  n. 

It  is  easy  to  compute  k^ (X).  However,  if  X is  unknown  and  one  applies 
tin-  Hayes  adaptive  procedure,  wi th  prior  gamma  «*(-^,\)  distribution  then, 

the  "desired"  ordering  level  at  the  beginning  of  tin-  i-th  month  is 

0 i-1 
K (T^)  - ^ where  T,  1 = ^ X-,  and 

i=0 
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-> 


(2.3) 


K^('l'i  -j)  “ least  non- negative  integer  k 


satisfyj ng:  G(kj2^. 


T._x)  - KG(kU.,v  + T._1)  * C 


in  which  G(k|i|i,v)  designates  the  c.d.f.  of  a negati\ e-binomial  distribution 
with  p.d. f. 


(2.1) 


0 < < 1,  and  \Aiere  i|r  = t/(1  + it).  The  "desired"  ordering  level 

be  always  attained  since  may  be  smaller  than  k\t^  ^).  We  thus 

(2.5)  Y?  = min(K°(Ti_1),Si)  - ^ , i = 1,1, ...  . 


cannot 

define 


The  upper- echelon  ordering  policy  is  defined  as  the  least  non-negative  integer 
Z such  that  He  predictive  probability  of  {(s^  + Z^  ^ - X^)+  - X^+x  + Z £ 
Ki+r,(Ti_i  + X.  + X.+x)),  given  T\  x,  is  at  least  y = 1-0.  We  have  estab- 
lished in  [1  ] that  the  Bayes  adaptive  ordering  policy  of  the  upper  echelon 
satisfying  the  above  requirement  is 


(2.6) 


0,  V 

Z.(T.  = least  non-negative  integer  Z satisfying: 

SSI 

G(Z  - -^(T^U^v  + T._x)  + Eg  (jH.,v  + T._1)[G  (Z-kJ(T1w1)  + 

+ Si  " Jl*i+l'v  + Ti+1  + " G (Zr  KiiTi-l)l*i+l'v  + Ti-1  + ‘1)]  * V 


wirere  t>^  — n + . > i - 1> . # • • • . 

In  Table  1 we  present  a simulated  two-echelon  system  with  parameters 
\ = 10,  v = 9*75*  t = 1.1,  c = .5[$],  p = 10. [ $],  c*  = 5. [$]  and  y - .85.  We 
denote  by  IQ,  IB  and  IX  the  values  of  and  X.. . KOP  designates 

K?(Tj_x),  IOR  designates  Y?  defined  in  ( C. 3 ) and  IZN  is  equal  to  Z*?^  ^) 
defined  in  (2.6).  The  cost  is  designated  by  CST  and  is  equal  to 
C*  • IZN  + C • max(0,  IQ,  - IX)  - p • min(0,  IQ  - ix).  Figure  1 shows 


- ( - 
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a (computer)  graph  of  the  moving- averager,  of  the  simulated  monthly  cost.,  CST. 

It  is  interesting  to  notice  that  in  the  present  example,  the  initial  stock 
level  and  S.  are  small  compared  to  the  "desired"  stock  level, 

which  is  k^(io)  = 28.  Although  X is  unknown  then  is  a relatively  quick 
convergence  of  K?('i\  ^ ) to  k^(10).  There  is  also  an  immediate  correction 
of  the  stock  level  f ^ by  a large  order  As  seen  in  Figure  1,  the  moving- 

averages  of  dST  converge  rapidly  to  the  limit  01'  their  expectations.  As  shown 
in  tlie  appendix,  if  the  system's  stock  level  g stabilises  on  a level  higher 
tlian  k" (X)  the  limiting  value  of  the  average  monthly  cost  is  approximately 


(2.7)  c = Xc*  + (c+p) [Q  2os(v;|X)  - X Pos(q-l|\)]  - p(Q-X)  , 

where  Q = k^(X)  - \.  For  the  parameter sof  .able  1,  q 13  and  c = 54.93. 

As  descussed  in  [ 4 ],  there  is  a high  correlation  in  the  simulated  data, 
between  the  order  level  Z.  and  the  demand  level  at  the  previous  month,  X.  , . 
Indeed,  in  Table  1 the  correlation  between  Z and  for  all  i a 3 is 

.')r>4.  it  was  therefore  recommended  in  [ 4 ] to  determine  Z^  and  Z„  exactly 
(according  to  (2.6))  and  to  set  / - x.  1 for  every  in?.  This  policy  will 
be  called  the  "approximate  upper-echelon  ordering  policy."  In  Table  2 and 
Figure  2 wo  present  the  simulated  system  for  the  same  demand  values,  X^, 
as  in  Table  1 v at  with  the  approximate  upper- echelon  policy.  We  see  that  the 
results  are  quite  similar  and  the  moving-averages  of  dFT  converge  to  the  same 
limit.  lie  convergence  is  however  somewhat  slower.  he  graph  in  Figure  2 is 
above  that  of  Figure  1 for  i £ 4o. 

?.  The  . ffects  of  the  Prior  Parameters. 

In  the  present  section  we  study  the  effect  of  the  prior  parameters  on 
the  control  procedure.  In  order  to  test  the  possible  effect  of  an  apparently 
wrong  choice  of  the  prior  parameter  we  have  performed  a simulation  run  similar 
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to  that  of  Tables  1 and  2,  but  with  the  prior  parameter  v 1.0  rather  than 
v = Jr.  '['able  3 wo  present  the  results  with  a::  exact  upper  -echelon 

policy,  wldle  in  Table  4 we  show  the  influence  of  the  approximate  upjvr 
echelon  policy.  As  seen  in  Tab_e  3 and  Figure  3,  the  inventory  system  adjusts 
itself  very  fast  to  the  difference  between  the  assunea  prior  expectation  of 
\ (which  is  l.l)  and  the  actual  value  of  \.  The  moving  averages  of  CBT 
converge  to  the  same  limit  value  as  those  of  Tables  1 and  : . Itiis  is  not  the 
case,  however,  when  we  apply  the  approximate  upper-echelon  policy,  as  seen 
in  Table  4 and  Figure  4.  One  has  to  apply  the  exact  upper- echelon  policy  for 
about  n - 10  months  before  switching  to  the  approximate  policy,  namely, 
using  for  T. . . Fxact  ordering  for  n = 2 mentns  only  is  insufficient 

due  to  the  effect  of  the  low  initial  stocks  and  relatively  high  demand.  As  a 
result,  the  Q.  values  in  Table  4 are  often  close  to  sero.  The  general  indi- 
cation is  that,  as  long  as  the  stock  levels  and  S are  too  small  one 

should  apply  the  exact  Bayes  ordering  policy  of  the  upper- echelon.  Once  the 

system  stabilises  with  S.  values  greater  than  K. (i.  .)  values,  one  can 

i i l-l 

switch  to  tiie  approximate  policy. 

4 . The  Effects  of  Non-Stationarity. 

Non- stat: onar ity  can  manifest  itself  in  different  forms.  We  investigate 
here  the  effects  of  sudden  unexpected  changes  in  the  intensity  of  the 

Poisson  process  of  demand.  More  specifically,  we  present  in  Table  ■s  a simu- 
lation of  an  inventory  system  with  = 10  for  Isis  10$  \ . = is  for 
"0  s i s 39  and  - 20  for  40  s i.  Ihe  statistician  is,  however,  unaware 
of  these  abrupt  changes  in  \ and  continues  to  use  the  Bayes  adaptive  policy, 
with  an  approximate  upper- echelon  ordering.  As  expected,  since  the  Bayes 
procedure  is  adaptive,  it  reacts  gradually  to  the  fact  that  the  observed 
demand  tends  to  be  larger  after  the  change  in  \ than  before. 


As  a consequence. 
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we  Liee  in  ’.’able  •>  that  the  K.(T.  , ) values  increase  steadily  alter  i = 20. 
They  do  not  reach,  however,  the  optimal  values  k^(l'>)  = 1*1  and  k^  (20)  bO. 
In  Figure  b we  see  the  graph  of  the  moving- averages  of  the  monthly  cost,  COT. 
We  cannot  compare  this  graph  with  that  of  Figure  2 when  i a 20.  This  is 
due  to  the  fact  that  the  expected  monthly  costs  associated  with  the  optimal 
policy,  basec  on  a complete  knowledge  of  the  epochs  and  magnitude  of  change 
in  X,  jumps  at  every  epoch  of  change  to  a new  steady  state  limit.  We  can 
compare,  however,  the  results  presented  in  Table  5 fnd  Figure  •'  to  those 
obtained  when  the  statistican  knows  the  initial  value  of  X and  is  not  aware 
of  possible  changes  in  X.  In  other  words,  the  same  optimal  policy,  based 
on  the  exact  initial  X,  is  employed  all  the  tine.  Ir.  Table  o and  Figure  6 
we  present  th  .•  simulation  results  of  such  a case.  Wo  see  in  Figure  6 that 
the  moving- averages  of  the  CST  are  below  those  of  Figure  b for  1 s i s E0 
(as  expected),  but  for  i > 20  the  adaptive  Bayes  procedure  performs  better. 

This  indicates  that  when  there  is  a possibility  of  ron- stationarity  a non- 
adaptive  optimal  policy  based  on  the  assumption  of  stationarity  is  likely  to 
perform  worse  than  an  adaptive  procedure,  which  corrects  itself  according  to 
the  observations  on  the  actual  demand.  We  remark  here  that,  as  seen  in  Table 
5,  the  Bayes  adaptive  policy  applied  reacts  rather  slowly  to  unanticipated 
abrupt  changes  in  X.  Indeed  this  policy  is  not  optimal  in  such  situations, 
since  it  has  been  designed  for  stationary  cases  in  which  X is  unknown  but 
fixed.  The  derivation  of  a Bayes  adaptive  policy  for  cases  of  abrupt  changes 
in  X,  at  unknown  epochs,  is  an  important  subject  for  future  research. 
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5*  The  Effect  of  Erroneous  Demand  Models. 

One  of  the  important  components  of  the  parameteric  control  model  is  the 
form  of  the  demand  distribution.  In  the  previous  sections  we  studied  various 
characteristics  of  control  procedures  derived  for  the  model  of  Poisson  dis- 
tribution of  demand.  It  is  interesting  to  investigate  how  robust  are  these 
adaptive  control  procedures  against  deviations  from  the  Poisson  model.  An 
exhaustive  suudy  of  the  sensitivity  of  the  Bayes  adaptive  procedures  requires 
more  development  in  terms  of  various  alternatives  to  the  Poisson  distribution 
than  what  we  present  here.  We  focus  attention  in  the  present  study  only  on 
the  following  alternative.  We  consider  a compound  Poisson  distribution  in 
which  the  mean,  X,  is  uniformly  distributed  over  the  interval  [10,  20]. 

This  provides  a discrete  distribution  with  a p.d. f. 

, 20  . . 

(5-1)  f(j)  = ^7JT  J*  e dX  = ^~[Pos(j|lO)  - Pos(  j [20)],  J =0,1,...  . 

00 

Notice  that  72  (l-Pos(j|X))  is  the  expected  value,  X,  of  the  Poisson  random 
0=0 

00 

variable.  Hence,  72  f ( j ) = 1.  Moreover,  the  expected  value  and 

3=0 

variance  of  a random  varialbe,  J,  having  the  above  p.d. f.  is  E{A)  = 15  and 

(5.2)  Var{ J}  = E[Var{ jJ  A) } + Var{E{j]A)}  = E(A)  + Var(A)  = 2.5  + ^ = 23.33.. 

Accordingly,  we  may  try  to  compare  the  behavior  of  the  Bayes  adaptive  control 
procedure,  for  a Poisson  distribution  under  a demand  distribution  specified 
by  (5-1)  to  that  of  a Poisson  with  mean  X = 15.  In  Table  7 and  Figure  7 
we  see  the  performance  of  the  Bayes  adaptive  control  under  the  compound 
Poisson  demand  distribution  and  the  exact  upper-echelon  ordering.  We  see 
that  Ki(Ti_1)  approaches  indeed  k°(l5)  = *+l.  The  upper- echelon  policy  of 
exact  Bayes  adaptive  control,  according  to  (2.6),  guarantees  that  the  system 
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stock  level,  if,  is  generally  above  the  "desirable"  level  K . ( T . ^ ) . The 
moving- averages  of  the  monthly  cost,  COT,  converge  consequently  to  a limiting 
expected  cost,  which  is  approximately  equal  to  100. ‘ . Notice  that  if  the 
demand  distribution  were  Poisson  with  mean  X - l‘>  then  the  limiting  expected 
cost,  according  to  (2.7)  would  be  7?  81.71*  Thus,  there  is  2 rj/&  increase  in 

the  limiting  expected  cost  due  to  the  fact  that  the  actual  demand  distribution 
is  not  Poisson.  It  is  of  interest  to  test  how  wsuld  the  approximate  upper- 
echelon  ordering  policy  perform  in  the  present  case,  compared  to  that  of  the 
exact  upper-echelon  policy.  In  Table  8 and  Figure  8 we  present  a simulation 
with  the  approximate  policy,  parallel  to  that  of  Table  7*  We  see  that  in  the 
case  of  an  approximate  upper- echelon  policy  tlie  system  stock  levels,  0 , 
are  generally  below  the  "desirable"  one,  aid  shortages  in  the  lower- echelon 
are  prevalent.  The  limiting  expected  monthly  cost  is  increased,  as  seen  in 
Figure  8 to  138.  This  is  about  3'$  increase  over  that  of  the  exact  upper- 
echelon  policy . Tlie  indication  is  that  tlie  approximate  upper-echelon  ordering 
policy  isnot  robust  against  considerable  deviations  of  the  actual  demand  dis- 
tribution from  the  assumed  one.  It  is  safer  to  use  .he  exact  policy- 


(n  Nummary  and  Conclusions. 

In  tlie  present  study  wo  have  shown  that 

1.  Tlie  exact  Bayes  adaptive  ordering  procedures  adjust  very  quickly 
to  the  actual  demand,  even  if  the  initial  prior  parameters  arc 
wrongly  chosen  or  if  the  initial  stock  values  are  inappropriate. 
The  moving- averages  of  the  actual  monthly  costs  rapidly  converge, 
in  the  stationary  ease,  to  tlie  limiting  ixpected  monthly  cost. 

2.  When  tlie  system  stock  level,  , is  above  tlie  "desired"  stock 

level  hb(T^  ^ ),  the  approximate  upper- echelon  ordering,  • X ^ 

provides  an  effective  simplification  without,  much  additional  cost. 
Op  the  other  hand,  if  the  system  stock  level,  S^,  is  too  low  the 
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upper-echelon  ordering  should  be  according  to  the  exact  Bayes 
predictive  policy. 

3.  In  cases  of  abrupt  changes  in  the  demand  distributions  at  unknown 
time  points,  it  is  better  to  apply  adaptive  procedures  tluui  optimal 
procedures  based  on  the  assumption  of  stationarity  of  demand. 

U.  'Hie  exact  Bayes  ordering  policy  protects  better  than  the  approxi- 
mate one  against  possible  discrepancy  between  the  model  ai.d  the 
actual  demand  distribution. 

Several  questions  discussed  previously  require  further  research.  In  particular, 
the  development  of  optimal  adaptive  policies  for  non- stationary  demand.  We 
have  applied  lore  the  Bayes  adaptive  procedures  designed  for  stationary  lV»isson 
demand  with  unknown  mean,  X.  If  X changes  at  unknown  time'  points  those 
policies  are  not  optimal  any  more.  How  much  more  complicated  the  optimal 
policies  would  be"  These  are  open  problems  for  further  research. 
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TABLE  4:  Simulated  rrwo-Echelon  System  under  Poisson  Demand, 
\ = 10,  Bayes  Control  v = 1.0,  t = 1.1,  y = .85; 
Approximate  Upper  Echelon  Ordering. 


I 

Iw. 

IS 

IX 

KOP 

JOB 

IZN 

CST 

1 

5 

1 0 

4 

7 

2 

Z> 

10.50 

Zf 

6 

if. 

1 1 

3 

12 

30.  0 0 

0 

Zf 

u_ 

14 

14 

170. 00 

4 

0 

0 

14 

20 

14 

210.  0 0 

5 

0 

0 

1 0 

24 

j 

14 

170. 0 0 

6 

0 

4 

10 

24 

4 

1 0 

150. 00 

7 

0 

3 

11 

25 

•I; 

1 0 

160.  00 

3 

0 

7 

15 

26 

1 1 

205. 00 

9 

0 

13 

27 

Zf 

15 

205. 00 

1 0 

0 

o 

6 

23 

o 

13 

125. 00 

1 1 

o 

3 

1 0 

27 

31 

O 

130. 00 

12 

0 

12 

12 

27 

12 

1 0 

170. 0 0 

13 

0 

6 

“7 

27 

6 

12 

1 3 0 . 0 0 

1 4 

0 

3 

1 0 

07 

3 

7 

135. 00 

15 

0 

1 1 

13 

07 

1 1 

1 0 

1 3 0 . 0 0 

16 

0 

5 

9 

5 

13 

1 55 . 0 0 

1? 

0 

5 

27 

9 

35 . 0 0 

13 

1 

14 

“7 

27 

13 

5 

35.  00 

13 

7 

16 

10 

26 

9 

7 

65.  00 

20 

z> 

11 

10 

26 

n- 

1 0 

30.  0 0 

21 

i 

3 

13 

26 

“7 

1 0 

170.  0 0 

■z>  z> 

0 

5 

1 0 

07 

3 

13 

165. 00 

23 

0 

■5 

14 

07 

3 

10 

130. 00 

24 

0 

4 

13 

27 

u 

14 

200. 00 

25 

0 

1 

13 

23 

1 

13 

135. 00 

26 

0 

3 

23 

Zf 

13 

145.0  0 

27 

0 

l 

14 

23 

7 

9 

130. 00 

23 

0 

5 

10 

23 

14 

170.0  0 

23 

0 

4 

6 

23 

4 

1 0 

110.0  0 

30 

0 

12 

3 

23 

12 

3 

110. 0 0 

31 

* 

14 

1 1 

27 

10 

3 

110.  0 0 

32 

9 

3 

23 

O 

1 1 

105. 00 

33 

1 

9 

1 0 

0- 7 

1—  1 

3 

3 

13  0.  0 0 

34 

o 

10 

10 

27 

1 0 

1 0 

150.  00 

35 

0 

3 

10 

27 

3 

10 

150.  00 

36 

0 

3 

10 

Zf' 7 

3 

10 

150. 00 

37 

0 

3 

*7 

07 

0 

1 0 

120. 00 

33 

1 

1 1 

1 0 

27 

1 0 

7 

125. 00 

33 

1 

1 1 

“7 

27 

1 0 

1 0 

110. 0 0 

40 

4 

1 1 

13 

27 

■7 

7 

1 

125. 00 

41 

0 

3 

11 

Zf~> 

u < 

3 

13 

175. 00 

42 

j 

4 

10 

0*7 
u-  1 

4 

1 1 

155. 00 

43 

fl 

-7 

1 

16 

07 

•7 

10 

210.  00 

44 

0 

Zf 

u 

3 

23 

16 

170.  00 

45 

0 

3 

3 

23 

3 

3 

75 . 0 0 

( "able  •*,  oont. ) 


l -AO** 


1a; 

IX 

KOI' 

IOK 

IZN 

.V  •' 

" 

— 

— 

4? 

;i 

t 3 

7 

3 — » 

4 ^ 

3 

::  5 . 0 0 

4 7 

A 

1 3 

J 

*!* 

35.5  0 

4 3 

1 0 

1 3 

3 

“* 

"\ 

3 

41 . 00 

4 3 

*5 

12 

13 

3 “!* 

3 

1 3 0 . 0 0 

5 0 

) 

1 

o 

1 

13 

155. 00 

51 

0 

3 

% 

3 -!• 

3 

3 

3 0 . 0 0 

5.3 

) 

13 

1 0 

3 “» 

1 3 

■> 

3 

I 15. 00 

53 

15 

1 0 

3*? 

■> 

1 0 

3 0.  0 0 

5 4 

*5 

3 

3 

3 "» 

■> 

3 

1 0 

j0.  0 0 

55 

> 

12 

1 0 

3 *5» 

t 0 

3 

no.  oo 

53 

> 

4 3 

4 3 

A 

3 — * 

1 0 

1 0 

15  0.  00 

57 

1 3 

3 

3 

12 

1 3 0 . 0 0 

53 

0 

3 

l l 

3 -!» 

3 

1 3 

1 '5. 00 

:•* I\«  1 : Mu' i.n*r,  Average  oV  Monthly  u\’ 

V 1,  T - 1.1,  p 101  $],  0 

:'xact  l^tpor  Koholon  0i\iering. 


v'OIlt  tv'l. 


C'*i.  00  #•.*. 


'.00  ♦> 


rV.OO  ► 


(‘>.00  *• 


*»♦♦♦♦♦♦♦♦♦ 


0.00  f * 


. . * * » * *.  . , 

t 0.  00  20.  00  30.  00  4 0.  00  50.  00 


r-;o4 


..  L : 


7.  imulatcd  IVo-eheloti  ;'.v.;tem  Uiu  . ol. 
Inoivtuun^,  Intensity;  10  (L  s l’i)/ 

X:  ■ ‘.'0  t i a-'<0)$  i'.'iyc;:  C 'ontrol  h 7 

Approximate  Upper  Xchelon  . ’rdor  li  if". . 


:-A>n  Demand.  With 
\i  - 1-  ('  0 s i s 
, T 1.1,  V . 


r 

In' 

17 

IX 

KOF 

12.:: 

— — 

« 

1 

*5 

1 0 

4 

3 1 

c 

J 

cr  “1 

26  0.  5 0 

3 

2* 

2* 

2» 

2 4 

J 

;i 

3.  00 

j 

0 

52 

1 4 

3w 

2 2 

2* 

170.  0 0 

4 

’3 

3 3 

1 4 

26 

13 

1 4 

1 3 0 . 0 0 

5 

* 3 

30 

1 

1 J 

w '3 

4 

J.  2* 

1 4 

71.0 ' J 

o 

3 

3 4 

1 0 

23 

1 0 

1 0 

5 4.  00 

y* 

13 

3 3 

1 1 

kl  7* 

1 0 

1 0 

5 3.5  0 

3 

i 

3 1 

15 

23 

1 1 

1 1 

5 6.  0 0 

3 

1 3 

22, 

1 3 

3 0 

1 7 

15 

"’5 . 0 0 

1 0 

1 •* 

3 0 

2* 

3 0 

1 3 

1 3 

7 0.5  0 

1 1 

2 4 

39 

1 0 

29 

rr 

j 

2* 

37.  0 0 

12 

13 

42 

12 

29 

1 0 

1 0 

5 3.5  0 

1 3 

1 7 

3 2* 

29 

1 2 

4 7‘ 

3*5.  0 0 

1 4 

“•  7* 

39 

1 0 

2 3 

2* 

—* 

41.0  0 

15 

1 3 

41 

1 3 

*2! '3 

1 0 

1 0 

52.5  0 

16 

15 

35 

2? 

1 4 

1 3 

2>  O*  . >J  'J 

17 

2 0 

3 2* 

5 

29 

•3* 

52. 5 0 

13 

2 4 

44 

2 3 

4 

5 

2*  3 • *5  0 

19 

*. 

*2  l 

46 

1 0 

7*  *7* 

4 0.5  0 

3 3 

1” 

41 

13 

1 0 

1 0 

3 0.  0 0 

21 

•3 

3 0 

15 

~t  -N 

LL  2* 

19 

13 

150.  0 0 

TO 

1 3 

25 

19 

29 

4 ~i 

I i* 

15 

1 35 . 0 0 

2 3 

2* 

24 

7»  7* 

30 

l ? 

19 

255. 00 

2 4 

3 

17 

1 4 

31 

4 C 

1 .* 

7*  7* 

2 3 0.  0 0 

~» er 
*2  J 

3 

7» 

4 y 
i • 

31 

19 

1 4 

210. 00 

26 

5 

7*  7* 

12 

3ii 

7*  7* 

17 

155. 00 

15 

29 

12 

2?  2 

1 4 

12 

31.50 

J a 

▲ 1" 

3 4 

15 

32 

15 

1 7* 

k *2 

6 1 . 0 0 

2? 

1 7 

31 

1 4 

-»  “1 

2**2. 

1 4 

15 

76.5  0 

3 0 

1 7 

29 

17 

32 

12 

1 4 

"7  |"|  ]j  ;"| 

31 

4 “< 

4 *_ 

7>  *r» 

16 

-«  -I 

2»  *2. 

15 

17 

1 25 . 0 0 

32 

1 1. 

25 

12 

3 3 

1 4 

16 

3 0.  0 0 

2*  2? 

1 3 

30 

13 

33 

17 

12 

1 1 0. 00 

3 4 

12 

23 

16 

33 

4 

1 2 

13 

130. 00 

3*5 

1 _L 

24 

23 

2*  3 

12 

i 3 

190.  00 

j 2 

1 

19 

4 “1 

1 il 

3 4 

12 

23 

225. 00 

37 

23 

t 4 

3 4 

16 

12 

130.  0 0 

33 

•3 

32 

1 1 

3 4 

2 2* 

1 4 

9 0.  0 0 

33 

21 

33 

7*  7* 

3 4 

12 

1 1 

35.  0 0 

40 

1 1 

25 

1 0 

34 

1 4 

110.50 

41 

15 

26 

2 0 

34 

1 1 

1 0 

1 00. 00 

42 

2 

23 

17 

35 

7*  7* 

7»  !*| 
kl  V 

210.  00 

43 

1 1 

21 

24 

35 

1 0 

17 

215. 00 

44 

A 

4 ~~y 

A * 

25 

35 

17 

2 4 

37  0.  00 

45 

0 

9 

23 

36 

4 

25 

355. 00 

J 


S'  'J1  .J» 


T-404 


-22- 


TABLE  6:  simulated  Two-Echelon  ;3ystem  Under  Eoi..on  Demand  With 

Increasing  Intensity;  X-j_  = 10  (i^l'y),  Xj_  = l'j  (20£i^39), 
= 20  (i^Xo);  Lower  Echelon  Control  Optimal  For  Poisson 
with  X = 10;  Approximate  Upper  Echelon  control. 


I 

IQ 

IS 

DC 

KOP 

1QR 

IZM 

CUT 

1 

rz 

_l 

« *. 

1 ■-* 

4 

3'  3j 

rr 

4 3 

215.5 0 

z> 

3* 

5, 

3*  '3* 

j 

j 

0.  0 0 

3 

0 

4 3 

1 4 

3*0 

23 

o 

17  0.  0 0 

4 

1 4 

29 

1 4 

23 

J 4 

14 

70.  0 0 

5 

1 4 

21 

10 

23 

7 

1 4 

72.  0 0 

Z‘ 

1 i 

25 

10 

3*0 

14 

1 0 

5 0 . 5 0 

7 

15 

29 

1 1 

3*0 

13 

10 

52.  0 0 

3 

4.  • " 

23 

15 

OO 

11 

1 1 

53.  0 0 

3* 

4 

1 J> 

23 

1 3 

23 

1 0 

15 

75.  00 

1 0 

10 

21 

j 

v3  -3* 

1 1 

13 

37.  00 

11 

15 

30 

1 0 

23 

13 

5 

32 . 5 0 

12 

13 

33 

12 

3**3j 

1 0 

1 0 

53.  00 

i 

4 J 

1 3 

3*  “7 

“7 

3*0 

u_ 

1 1 

1 O 

4 i— 

34.50 

1 4 

20 

30 

10 

23 

o 

7 

4 0.  0 0 

15 

13 

3|'3* 

13 

w'3* 

1 0 

1 0 

52 . 5 0 

IS 

15 

23 

•j 

1 1 

13 

S3.  00 

1 “7 

4 1 

17 

3*7 

5 

23 

1 0 

•3 

51 . 0 0 

13 

3'  3* 

. — L. 

35 

23 

5 

1 j :j 

13 

21 

~;"7 

10 

23 

7 

“7 

4 0 . 5 0 

20 

13 

32 

13 

23 

1 0 

10 

5 0.  0 0 

21 

1 0 

21 

15 

23 

1 1 

13 

14  0.  0 0 

3*2 

“i 

1 Z> 

13 

3*Q 

1 0 

15 

2 05.  0 0 

“* 

0 

15 

3*  3* 

23 

15 

19 

315.  0 0 

24 

0 

3 

14 

23 

j 

"•  Z' 

250. 00 

25 

0 

13 

17 

yy 

13 

14 

240. 00 

26 

0 

13 

12 

23 

13 

17 

205. 00 

27 

i 

20 

12 

23 

1 4 

12 

12  0.  00 

3*  Q 

>—  'j 

ZJ 

25 

15 

23 

17 

1 “• 
i iZ 

13  0.  0 0 

23 

1 0 

3*  3* 

14 

3*  ;3* 

12 

15 

115. 00 

30 

s 

20 

17 

3*  O 

12 

1 4 

130. 00 

31 

3 

13 

13 

3*0 

15 

17 

215. 00 

3*  3* 

Ou 

13 

12 

23 

14 

IS 

1 3 0 . 0 0 

33 

4 

21 

13 

23 

17 

12 

2 0 0.  00 

34 

13 

13 

23 

13 

13 

220. 00 

35 

3 

15 

23 

23 

12 

1 '5 

23  0.  00 

33 

0 

10 

12 

23 

10 

23 

235. 00 

37 

0 

14 

14 

23 

14 

12 

200. 00 

33 

•j 

23 

11 

23 

23 

14 

130.  00 

33 

12 

24 

22 

23 

12 

1 1 

155. 00 

40 

o 

13 

10 

23 

1 4 

~J  O 

L—  L_ 

130. 00 

41 

.3 

17 

20 

o o 

1 1 

10 

13  0.  00 

42 

0 

13 

17 

23 

13 

20 

270. 00 

43 

12 

24 

23 

10 

17 

305. 00 

44 

0 

3 

25 

23 

'-! 

24 

370. 00 

45 

0 

0 

23 

23 

0 

25 

355. 00 

I 

IQ 

IS 

IX 

KOP 

ICR 

IZN 

CST 

J 

43 

0 

1 

25 

23 

1 

23 

335. 00 

, 

47 

0 

1 

19 

23 

1 

25 

315. 0 0 

43 

0 

cr 

_i 

17 

23 

5 

19 

235. 00 

43 

0 

13 

17 

23 

l - 

17 

255. 00 

, 

50 

0 

15 

17 

23 

15 

17 

255. 00 

\ 

51 

0 

15 

13 

23 

15 

17 

275. 00 

j 

52 

J 

13 

19 

23 

13 

19 

235. 00 

• 1 

53 

11 

23 

OC; 

1 1 

19 

325. 00 

■' 

54 

0 

”7 

31 

Q 

■7 

23 

425. 00 

55 

0 

0 

14 

23 

0 

31 

295. 00 

53 

0 

9 

-O 

t—  L- 

23 

9 

14 

290. 00 

57 

0 

13 

23 

23 

13 

22 

34  0.  00 

53 

0 

9 

13 

23 

9 

23 

295. 00 

; 1 

FIGURE  6:  Moving  Averages  of  Monthly  Cost;  Fixed  Control 
Under  Poisson  Demand  With  Increasing  Intensity. 


28‘  .00  <■.  ♦ » *. *... ♦.  , * 

-I 

-I 

-I 

22^.00  > +I 

’ ♦ -I 

-I 

*♦♦♦♦  -I 
♦♦♦♦♦♦♦  - 1 

150.00  f 

* *♦♦♦♦♦♦  _I 

* - 1 

* * ♦♦♦♦♦♦♦♦  _j 

•■♦♦♦♦  _ 1 

75-00  ► >♦  fj 

»♦  -I 

-I 

-I 

0.00  * » 'I 

• * ♦ ♦ * » », . 

0.  00  10.  00  20.  00  30.  00  40.  00  50.00  30.00 


i 404 


TABLK  Y 

: I'Lmul 

abed  IVi'-  vln 

•Ion  :’>y tern;’ 

; Under  0 

otnpound  1 

i 'oi.:*.  :>OM 

Pemmul  With  Kandon 

i Un.il\)rin  [ntennity 

( 10  < x 

5 '0); 

Bay  on 

Control  v 

'’■Y1',  T 

1.1,  V 

• i 

Kxaet. 

Upper 

ivhelon  Ordering. 

1 

B! 

IB 

IX 

KOI’ 

I0K 

i.'X 

' • * .1 

" 

' 

1 

5 

1 0 

2 

7 

2 

2 

2 0.  0 0 

*> 

1 

4 

1 3 

15 

2 

20 

25  0.  0 0 

j 

o 

0 

1:3 

2 4 

0 

25 

305. 0 0 

\ 

o 

2 

1 3 

30 

2 

3 0 

230. 00 

*5 

o 

1 4 

12 

31 

1 4 

15 

195. 00 

2 

15 

31 

29 

12 

1 9 0 . 0 0 

13 

32 

1 3 

3 3 

13 

13 

91.5  0 

2 

19 

31 

13 

3 3 

12 

1 3 

^5 . 5 0 

1 3 

31 

13 

3 4 

13 

21 

1 35. 00 

i o 

15 

23 

1 3 

35 

1 3 

13 

31.  0 0 

1 1 

15 

33 

17 

35 

20 

1 3 

35.  0 0 

12 

13 

^ ,■ 

2* 

3 o 

13 

13 

93 . 0 0 

1 3 

30 

4 4 

1 0 

3 4 

4 

2 

25.  0 0 

1 4 

24 

52 

21 

3 4 

1 0 

1 0 

51.5  0 

15 

1 3 

3 4 

13 

35 

21 

J 2 

1 35 . 0 0 

13 

13 

2o 

2 1 

“t  * 

2 2 

1 0 

19 

1 45. 00 

l 7 

*5 

23 

2 2 

U l_ 

3 1 ' 

23 

*>  ~i 

lL  9 

235. 00 

13 

o 

25 

12 

33 

19 

2 4 

1 3 0 . 0 0 

19 

1 3 

3o 

1 4 

33 

•s 

^ 2 

12 

7 0 . 0 0 

20 

Zf  Zf 

43 

1 1 

3 3 

13 

l 4 

75 . 5 0 

21 

J 

4 7 

12 

3 . 

1 0 

1 0 

57. 50 

2 2 

25 

49 

17 

3 , ■ 

12 

12 

34.  00 

2 3 

J !J 

42 

19 

3 »■ 

1 

i i 

1 7 

35 . 5 0 

2 4 

13 

-» rr 

2 j 

15 

33 

17 

2 0 

101. 5 0 

*>rr 

w _l 

2 0 

37 

1 4 

3 3 

17 

15 

73 . 00 

23 

ki  J 

4 3 

1 1 

33 

15 

1 4 

. o . 0 0 

O 

27 

47 

21 

33 

1 1 

1 1 

53.  0 0 

23 

17 

4 0 

25 

33 

21 

2 2 

1 0 • i . : 1 0 

29 

1 3 

23 

1 0 

39 

1 3 

i , - 

kl  *2 

1 31.50 

3 0 

13 

3 3 

1 3 

39 

*2  2 

k— 

1 0 

5 1 . 5 0 

31 

25 

51 

1 3 

33 

1 3 

12 

. 0 0 

32 

25 

43 

2 0 

33 

1 3 

1 3 

» 7 . 5 ' i 

3 3 

13 

4 0 

39 

21 

21 

1 10.50 

3 4 

32 

43 

17 

33 

2 

2 

37.5  0 

35 

21 

5 0 

1 1 

33 

17 

1 r* 
i . 

90.  00 

3 3 

2 

u.  i 

45 

23 

33 

1 1 

1 1 

57.  00 

2 

15 

39 

17 

39 

24 

24 

1 40. 00 

33 

:•  2 

t— 

j 3 

15 

39 

l 1 

17 

33.  50 

39 

13 

42 

2 ■? 

U 1 

39 

21 

15 

1 .->5 . 0 0 

$0 

12 

32 

13 

39 

20 

23 

13  0. 0 0 

41 

13 

31 

12 

39 

15 

13 

32.  0 0 

♦2 

19 

4,’ 

1 3 

39 

20 

12 

3.00 

♦ 3 

23 

5 0 

39 

1 3 

1 3 

"4.50 

4 4 

32 

55 

1 0 

39 

2 

41.0  0 

45 

29 

5 3 

13 

39 

1 0 

1 '1 

55.5  0 

(Table  7,  eont. ) 


13 

13 

J5. 

0 0 

1 1 

1 1 

j4. 

:.l  t.l 

"i 

•4 

5 3. 

1 !j 

1 3 

i 3 

7 3. 

S !.l 

4 0. 

0 0 

2 0 

21 

1 75. 

4 0 

1 \ 

23 

131. 

5:.i 

21 

43. 

0 0 

1 1 

1 1 

.->2 . 

0 0 

1 3 

1 3 

y?. 

0 0 

17 

1 7 

37. 

0 0 

17 

17 

91 . 

50 

3 

3 

44. 

5 0 

KIOllKK  7:  Miv'Lii(',  Awrap.ra  ol’  Monthly  dost.;  dotirjxnmd  Ibissoti  Pomand; 

Bayes  Oontro,  v 9.7‘',  t 1.1,  \ . 8';,  p 10[$], 

e •''l-t'],  c*  - ‘'|  $];  Kxnct  Upper  Kohelon  Ordering. 


oo  t.  *. 


oo 


10.00  f 


T’.oo  „ 


0 . iK)  * 


»»* 

* • 


* ♦ * * »♦♦♦«> 


* *.1-1 

-I 

-I 

-I 

M 

-I 

-I 

-I 

-1 

M 

-I 

► -I 

**♦♦♦♦♦♦♦♦♦♦♦♦  -i 

-I 

M 

-I 

-I 

-1 

-1 

•■1 


0.00  10.00 


, . . * »,  . , 

20.  00  30.  00 


♦ 0 . 0 0 50.  0 0 > 0 . 0 0 


. • t mu  1 at* \i  V*  l.e-rn.  Uiuit-r  ‘onpound  oi;:roti 

l '•■in;  it  ul  i t ) i i 1 1*  'in  ! • . L t\ ■<  rf.  Inti'ii.  i tv  ilOS\S'i))j 

■ AV'V'  v , i 1.1,  > . ; A]>vr>'xi- 

m: i I..  i I'm ti  • i*  ."I i 1 ■ ■ 1 1 i\l. -r  i lie.. 


I 0 
1 1 
12 


l 3 

i » 


i 3 


13 

13 

2 0 


21 


2 V 


-•'3 


23 

"V| 


l 1 


l 0 


■4 

1 J 


20 


1 4 


13 


13 
l 0 


0 


1;’ 

1\ 

KOI 

17..N 

1 J 

3 

2 7 1 • . i 

4 

1 3 

“•  “7 

0 

j 

1 2 . 0 

*0 

1 3 

32 

J*  w 

1 3 

2 3 0.  0 

l 3 

3 3 

1 3 

> i.5 

21 5 

12 

3 3 

1 

1 3 

35 . 0 

31 

15 

35 

1 0 

1 J 

3 i . r 

23 

1 3 

j ^ 

t *5 

15 

7 ? . 5 ' 

1 3 

_•’»  j 

12 

1 3 

35.  0 

25 

l 3 

3 

i r 

. J 

13 

150.  0 

22 

1 3 

J 7 

» 

1 > 

12  0.  0 

27 

1 *7 

i 

3 , 

1 3 

1 3 

1 45.  0 

2» 

2*  3 

1 3 

17 

37 . 0 ' 

3 3 

1 0 

t 3 

JV 

35 . J 

10 

21 

35 

12 

i :i 

5 1.0: 

-icr 
^ j 

13 

57 

' 1 

21 

1 45.  0 

l 7 

->  j 

3 , 

1 0 

13 

2 3 i . 0 C 

17 

1*  3 

3*3 

17 

■>  i 

w 1 

325. 0 • 

1 3 

1 2 

39 

1 3 

»2L  2 

2 3 0 . 0 

“»  j 

1 4 

39 

21 

12 

13  0.  0 

3 0 

1 1 

39 

1 4 

l 0 0 . 0 

31 

12 

33 

12 

1 1 

5 3.5'. 

3 3 

i 

L • 

33 

1 4 

12 

1 . 0 . 

2 1 

13 

3 3 

1 1 

17 

115.0 

20 

15 

39 

1 3 

13 

135. 0C 

*)  7« 

1 ♦ 

3? 

17 

15 

135.  0 

-*  “? 

1 1 

39 

1 ? 

1 4 

1 0 0.0 

31 

21 

33 

15 

1 1 

1 95 . 0 5 

24 

25 

3? 

1 4 

21 

255.  O' 

1 0 

1 0 

4 0 

1 0 

25 

225. 0 v 

21 

1 3 

39 

21 

l 0 

1 '3  0 . 0 ' 

33 

1 3 

39 

25 

1 3 

115.0: 

30 

2 0 

39 

1 0 

1 3 

.•>5.  J . 

23 

"7 

39 

1 3 

2 0 

i oi.5 : 

23 

17 

39 

1 3 

"7 

45.  0 

32 

l 1 

39 

L.  J 

t 7 

35 . 5 

23 

23 

39 

•«1 

1 1 

75.  0 

■?  *3 

w U 

1 7 

39 

1 7 

j* 

2 35.  0 0 

1 3 

15 

39 

1 1 

1 7 

135. o : 

2 1 

■>  *7 

39 

15 

3 35.  i 

l 1 

1 3 

4 0 

1 4 

■>7 

235. 00 

1 3 

12 

♦ 0 

1 3 

1 o 

200.  0 0 

23 

1 3 

4 0 

*7 

12 

1 3 0 . 0 •.» 

31 

"7 

♦ 0 

1 

1 3 

3 . 0 • 

3 3 

1 0 

39 

12 

7 

42.  0 i 

33 

13 

39 

l 3 

1 0 

5 4 . j i 

21 


1 


1 1 o . :•  -i 


31 


' UU'KO  fi;  Moving  Avorogt*! 

Payee  Control, 


oi 

v 


Monthly  Ooetj  'owjxuuul  i‘oi.:::on  emandj 


1.1, 


P 10|  ■ ], 


>|  41  ]}  Approximate  Op’x<r  .‘•'elu-lo.’i  Oial  fin. 


: \'v  >.00  1 1- . * 


* 


* 


* 


» 


* 


>.00  *. 


l‘)0.00  >. 


* 

» * ♦ • 


* * * »***♦ 
**♦*♦♦*♦ 


7' <.oo  *• 


o.oo  *.  » 

. . ♦ * * * *■ • •. 

3.  00  1 0 . 0 0 2 0 . 0 0 3 0 . 0 0 10.  0 0 5 0 .00  i 0 . 0 


T-404 


8.  References. 

[1]  Kaiffa,  H.  and  Behlaifer,  R.  (idol) 

Applied  statist. Leal  Decision  Theory.  Cambridge:  Harvard  University  pres; 


Packs,  B.  (ldoO ) 

Bayes  Sequential  Design  ol’  Stock  Levels,  Naval  Res.  Legist,  Ouart. . 
lo:  lU3-l‘)'p.  

Packs,  S.  (loyo) 

A Two- Kehelon  Multi-Station  Inventory  Model  for  Navy  Applications; 

Naval  Res.  Legist.  Quart..,  17 : 70-8‘>. 

Packs,  S.  and  Fennell,  J.  (1072) 

Bayes  Adaptive  Control  of  Two- Kehelon  Inventory  Systems,  I:  Development 
ioi  a special  ease  01  One-station  Lower  Kehelon  and  Monte  Carlo  Kvnluation; 
Naval  Res.  Legist.  Quart..,  10;  l'->-28. 

Packs,  S.  and  Fennell,  J.  (W73) 

Distribution  ol  Adjusted  stock  Levels  Under  Statistical  Adaptive  Control 
fl-ocedur :s  of  Inventory  Systems;  Jem-.  Amer,  Statist.  Assoc..  68:  88-01. 

Packs,  S.  and  Fennell,  J.  (197I+) 

Bayes  Adaptive  Control  of  Two- Kehelon  Inventory  Systems,  II:  The  Multi- 
station Case;  Naval  Res.  legist.  Quart.. , 21:  57  >-‘^3. 

sacks,  S.  (107)1 ) 

On  the  Optimality  of  the  Bayes  Prediction  Policy  in  Two- Echelon  Multi- 
station Inventory  Models;  Naval  Res..  Legist.  Quart. ; 21:  pod- ‘17)4 . 

Packs,  S.  (1976) 

Review  o statistical  I'roblems  and  Methods  in  logistics  Research,  Oh.  10 
in  Modern  Trends  in  Logistics  Research  (w.  H.  Marlow,  ed. ),  The  MIT  Press, 
Cambridge,  Mass. 


-29- 


T-404 


9-  Appendix. 

9.1  Derivation  of  the  Limiting  Average  Cost. 

The  expected  monthly  cost  associated  with  Y?  and  P?  is 

(9.1)  E(OST)  . C*  I . E|c[i,;.  ♦ y“  - >9)*  - 

-p[<V  yJ  - xt)+  - ii(1ri  . 

If  S.  i K.  (T.  ., ) then  Q.  + Y?  = K.(T.  .,).  Furthermore,  we  have  seen  that 
i i—l  i i 1 i—1  7 

after  the  system  adjusts  itself,  Z?  « ^ and  K (T_.  «;  k^(X).  Accord- 

ingly, the  limiting  expected  monthly  cost  is  approximately 

(9.2)  M(\)  ■-  c\  + E{c[(k°(X)-  X.)*  - X.+1]+  - p[(k°(X)  - X.)+  - X^f)  • 

If  we  replace  kLVx)  - X^  in  (9.2)  by  its  expectation,  Q,  = k°(X)  - X,  we 
obtain  formula  (2.6)  namely 


(9.3)  c = X(c*  + p)  + (c+p)[Q  nos(QjX)  - X Pos[ (v}-l|X) ] - p;J  . 

To  evaluate  this  approximation  we  develop  the  formula  of  M(X).  From  (0.2)  we 
obtain 


kv  k°  ’ 

(9.M  M(X)  = c"X  + c D p ( J 1 X ) D p(i|x)(k°  - j - i)  + 

i--0 


0 

LmJ 

j=0 


0 


+ P £ P ( J lx)  D p(i|x)(i-k°+j)  + Xp 


00 

V' 


j=0 


i=k°- j+1 


j=k°+l 


p(j|x)  , 


0 c 

where  k = k (X).  Notice  that 


(9.5) 

and 

(9.6) 


c 

V 


i=l 

os 


p(ijx)vc-i)  = c Pos(c|X)  - X Pos(c-ljx)  ■=  (c-X)  Dos(e-l|x)  + cp(c|x). 


(i-e)  p( i | X)  = (X-c)  - (X-c)  Pos(c-l|x)  v ep(c|x)  . 

i=c+l 


-JO- 


T-404 


Hence,  (9.**)  obtains  the  form 
(9.7) 


K°_l 
v> 


M(\)  = (c*  f p)X  ^ (c*p)  S p(j]\)  [(k°-j-X) 

j=0 


Pot*(k°-J-l|X)  f (k°-j)  p(k°-j|x)]  - 


p[(ku-X)  lX>ii(k'‘-l|X)  + k°p(k°|X)]  . 

Finally,  we  notice  that 


k°-l 

(9.8)  D p( J 1 X)  Pos(k  - j-l|X)  = P{X  * X0  s k°  - 1)  = Pos(k°  - l|2X)  , 

j=0 


where  X^  and  X„  are  independent  identically  distributed  Poisson  r.v. *c 
with  mean  X.  Similarly, 


k°-l 


k°-: 


(9.9)  Jp(j|.\)  Pos(kL-j-l|X)  = X L'  p(o(x)  Pos(k°-2-jfX)  = XPos(k°-: 

J=o  j=0 


|2X). 


Accordingly,  from  (0.8)  and  ('>.9)  we  prove  that 


(0.10)  M(X)  = X(cOp)  + (c+p)[(k°-2X)  Pos(k0-rjrx)  + (k°-X)  p(k°-l|2X)J  + 
+ (c+p)[k°p(k°|2X)  - k°p(o|x)  p(k°[x)  - Xp(k°-llcx)  + 

+ Xp(olx)  p(k°-l| X) ] - p(k°-x)  P08(kC-l|x)  - pk°  p(k°lx)  . 


This  function  is  considerably  more  complicated  than  c given  by  (2.0)  or 
V' 1 . o ) . In  the  following  table  we  give  a few  numerical  comparisons  between 
c and  M(X).  We  see  that  for  X not  exceeding  10 


Table  9.2: 


Numerical  Comparison  of  c and  M(X)  for 
X = )(5)30  and  c*  - 9,  c = .9  and  p = 10. 


X 

c 

M(X) 

5. 

28. 471*48 

28.0891 ' 

10. 

*s4. 92688 

■Jl.l4090 

19. 

81. 70487 

79.75215 

CO 

107.  '*0707 

109.  09482 

25 

133. 20884 

131. 1*0261 

30 

159.00808 

1‘X'.9U419 

J 


the  approximation  to  M(\)  given  by  c is 
greater  than  10  the  approximation  is  not 
for  X s 20  if  less  than  2$.  since  (9.3)  is 


very  good.  For  values  of 
as  good,  but  the  relative 
considerably  simpler  than 


X 

error 

(9.10) 


it  seems  justifiable  in  many  circumstances. 
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